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Abstract 

We propose a measurement on the elastic vrvr scattering phase shift 
difference 6q — 6q through — > J/^ vr+vr" process in future high statis- 
tics BES-III experiment. The decay amphtude is constructed with seven 
Lorentz invariant form-factors and is compared with their theoretical es- 
timation. It is found that the phase shift difference can be obtained, 
based on a Monte Carlo study and it is expected the phase shift in the 
energy region between 350 MeV to 550 MeV can be measured at future 
BES-III. 
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1 Introduction 

In recent years, the operation of a number of high precision, high statistics 
experimental machines, varying from fixed target experiments to B factories, 
opens a new era for precision hadronic experiments. Based on that, both ex- 
perimental and theoretical studies on low energy vrvr and irK system produced 
in production processes have also received revived interests. The importance 
of these studies follows from the fact that when final state theorem applies, 
one can extract low energy elastic scattering phase shifts in the related chan- 
nels through a partial wave analysis. The information on vrvr, vrA' phase shifts 
then provides a crucial ingredient in understanding the dynamics of goldstone 
bosons and the spontaneous breaking of chiral symmetry. 



The experimental and theoretical activities in the last few years mainly 
focused on semi-leptonic and hadronic D decays (see for example, Refs. [1] - 
[1]). It is known that in D semi-leptonic decays p-wave dominates, and the 
more interesting s-wave component is small. In this paper we re-investigate 
the TTTT final state interactions in the — > J/'^tt^'k~ process. Here s-wave 
dominates and the next contribution comes from the tiny rf-wave. The exis- 
tence of the latter is however crucial for exploring the s-wave phase motion 
through interference effect. The decay product under concern is a three body 
final state, the J /'^ particle is however irrelevant to any final state interactions 
here. Using color transparency argument it is understood that the effect from 
rescattering between the J/\l/ and one of the pions is negligible. Another im- 
portant fact is that, in the kinematic region under concern, between the initial 

and the final J/'^titt there is no other on-shell intermediate hadronic state 
available (or are double OZI suppressed and hence negligible). Hence the final 
state theorem is applicable to the tttt system in J/\E^7r+7r~ process. 

The \Ef' — >■ J /^>TT^TT~ process has been the subject of a number of previous 
pubhcations ([5]- [7]). In Ref. [8], the author proposed a method to extract tttt 
phase shift from J/\E'7r~'"7r~, similar to Pais-Treiman [9] method for ob- 

taining TTTT phase shifts from Ki^^ decays, but only considered three partial wave 
amplitudes for reducing the difficulty of the analysis. A similar method [10] 
was also proposed in T(3S') T{lS)Tm process but only the lowest order 
in the pion momentum expansion was considered. In this work we are able 
to provide a more general parametrization to the decay amplitude comparing 
with what is given in most of previous papers. Our parametrization will be 
discussed in section 2. Furthermore we will also provide a Monte Carlo study 
in section 3 to test the stability and reliability to use our parametrization to 
extract the phase shift data. 

2 General structure of the J/^tt+tt" de- 

cay amplitude 

2.1 The Lorentz invariant form- factors 

There are 3 independent momenta Pn+,Pn- ,Pj/^ = Ps, which can be re- 
expressed in 3 variables q = + PTr-,P = Pn+ — Pn-, and p^. The three 
independent momenta can form 2 independent Lorentz invariant products, 
chosen as and p ■ p^ here. Then, 




g - p = 0, 



(1) 
(2) 
(3) 
(4) 



q-p; = -{Ml,-Ml 




(5) 
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Figure 1: ^3 is the angle between the beam and J/\I/ in the rest frame. ^* 
is the angle between the J/\l/ direction and the vr"*" in the di-pion rest frame, 
is the azimuthal angle between the beam J/\l/ plane and tt+tt^ plane in the 
rest frame (not drawn in the figure). 



with the kinematics factor of dipion system p = \Jl — the energy and 
the momenta of J/\l/ in the lab frame (the rest frame), which £^3 = (M^, + 



M| — s)/(2M^/) and \ps\ = y — M| are functions of g^. Moreover, p-ps can 
be expressed in the rest frame with the variables in the dipion rest frame, 

P-P3 = 27ip;;i cos^:(/3E3 + m), (6) 

where \p*^\ = p^fsjl is the three momenta of vr in the dipion rest frame 
and 6'* is the angle between the direction of 7r+ and direction opposite to the 
final J/\l/ in the dipion rest frame, see Fig. 1 for illustration. (5 = ^ 



is the boost factor from the dipion rest frame to the lab frame with 7 = 
(M|, -M| + s)/(2M*,v^). 

Denoting the polarization vector of and J/\l/ by e' and e respectively, 
we can form five invariants bilinear in e'e*: 

(e-e')> (e-g)(e'-g), (e-p)(e'-p), (e-g)(e'-p), (e-p)(e'-g). 

Hence, the amplitude have the independent structure: 

T = (e.e')Fo + (e-g)(e'-g)Fi + (e-p)(e'-p)F2 + (e-g)(e'-p)F3 + (e-p)(e'-g)F4. (7) 

Fi are the functions of s and p ■ ps and these allow us to do partial wave 
decomposition according to 6**. 

2.2 Partial wave decomposition 

From Ref . [11] , we obtain the basis of tensors 

= 1, (8) 
i^'^ = P', (9) 

p) = p^f - 1^2^/^-, with = 9^" - (10) 
3 q'^ 



where every tensor i^^^ transforms irreducibly as a tensor of spin L. In the 
present problem we have the four vectors, q'^, P3, e''^ and e'^ , which are inde- 
pendent on 6*. We can build the independent Lorentz scalars: 

g^, q ■ e, q ■ e' and e ■ e'. 

The available Lorentz vectors would be 

„M ^'M 

With these we go first to build L = quantities. This can only be obtained 
through t'-''-* and taking into account that the polarization vectors e' and e must 
always be contracted at the end, 

S = . {Jf )(s)(e' ■ 6) + Jf (s)(6' ■ g)(6 ■ q)}, (11) 

which can be expressed through S = e'f^e'^S^^, with S^^, = h{s)g^u+his)q^qu/ s. 

The (i-wave is more complicated, since we have to use i'^^) the number 
of contractions gets larger. The available four-vectors, q'^,p'^,e''^,e^, are then 
contracted with t'^^-^ in all the different possible ways: 

D = t{5 ■ {l!t\s)ie' ■ e)q^q^ + jf (.)(e' ■ e)p',p^, + jf (.)(6' ■ e^Q^ + q^p^,) 
+/f (s)(e' ■ g)(6 ■ q)q^q'' + jf (.)(6' ■ g)(6 ■ q)p^,p^, + (s)(6' ■ g)(e ■ + g>3^) 

+4°)(.)(e ■ g)(6'^g'^ + g^e^) + /;?(s)(e ■ q){e"^p^, + p^^e^) 
+/ff)(.)((6' ■ g)(6V + g^e-^)) + /l2H^)(e' ■ q){e^p's + P'sen 

+li?is){e"^e'' + e'^e'n} . (12) 

It is not difficult to find that there are not actually so many independent 
Lorentz structure. It can be written in a more compact way through D = 
e'f'e'^D'^^, with 

D'^^ = /^(,)J-^M.[(p3 .p)2 _ \p'(p,Y] + /:(s)igV[(P3 -P)' - \p\Pl)] 

S O S O 

+/^(.)4bV(P3 ■ P) - \pVp'3] + I'e{s)^[p'q''{P3 ■ P) - V^39l 

S o S (j 

+r^[s)-yp'' - ^-p'r"]- (13) 

In order to avoid that any form-factor becomes artificially large os small, 
we extract Lorentz structures that are numerically order one. The amplitudes 
are then expressed in the form 

5^'^ = /i(s)(7'^- + J2(s)W, (14) 

D'^^ = /3(.)^?^^[C0S2 - ^] + /4(5)^gV[cOs2 e: - 1] 

+h{s)^j^yp''{p, ■ p) - \pVpI] + ^6(^)^574^b'^9^(P3 ■ P) - \p'pU^ 

+J7(s)-bV- V^n, (15) 

S 6 



Until this point the derivation is completely general. Now, we make the 
main assumption: we will assume that no further rescattering occurs between 
the J/\Ef and the vrvr system. Hence, the phase-shift of the amplitude is due 
to the TTvr final state interaction. This allows to use the Watson theorem for 
the elastic scattering region (from the practical point of view, up to the KK 
threshold). The decay amplitude can be then decomposed into partial- waves 
(5*, D...) with their phase-shifts equal to those in tttt scattering (respectively, 
(5o, 52...): 

T = e'^e''[S^,e'''+D^,e''% (16) 

The D-wave is supposed to be suppressed with respect to the J = component 
and higher partial waves are neglected. 



2.3 Theoretical estimation to the leading contributions 

Starting from the effective lagrangian in Ref. [I2], which is constructed using 
chiral symmetry and heavy quark symmetry, the amplitude has the form, 

+92[p:.p:-+p:-p:.Ke,}, (n) 

where = (i?^±,p^±) in the rest frame. We can calculate these leading 
terms' contribution to our form factors /^(s), 

h(s) = -|j[f (. - 2™^) + iifd - ^f)^,,i£ifl (18) 



2s ^ p{s) 
1 

1 



his) = -^g2il - ^), (19) 



his) = —g,p{s)'s(3'^', (20) 
his) = ^92 , (21) 



with the rest being vanishing at leading order. This calculation suggests that 
the form factors I^is), his), his) are small quantities and the theoretical 
prediction can be checked by future experiments. Recall that the fit to the 
data can also in principle determine his). 



2.4 Expressions for angular distribution 

For three body decays [13] , 

dr = 77^,T7^2\-^\'\p:^\\p^^^^ (22) 
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where (p^*^, ^l-^) is the momentum of 7r+ in the dipion rest frame, and Q j/^ is 
the angle of J/^ in the rest frame. Hence, we can express in the form: 

— — — — OC Gklm is) COS'^ 0; cos' 63 COS™ O*^ 

dsa COS old COS d-idd)* ,, 

TT o Tit klm 

+ COS C sin ^3 cos ^3 sin 91 cos (Go(s) + ^2(5) cos^ 91), (23) 

where k,l = 0, 2, and m = 0, 2, 4. Gkim{s) are functions of /j(s) and cos(5o— ^2)- 
If we can determine Gkim{s) experimentally, then we can determine Ii{s) and 
extract the phase shift difference 5q — 

Alternatively, we can find the information in the partial distribution. What 
we find in the partial distributions is of the form: 

d'^T 



dsd cos 9* 

d^T 
dsd cos ^3 

d'^T 



Ao(s) + A2{s) cos' 91 + ^4(5) cos^ 91, (24) 

5o(s) + 52(s)cos2^3, (25) 
Co(s) + C2(s)cos2 0, (26) 



dsd cos (f) 

and another weighted distribution, 

W[s, cos 91] = dct>l d cos ^3 cos </.; cos ^3 

J— TT J — 1 UoU COS C/^tt LOS (73Cl<)t'jj- 

y d COS 6^3 cos^ 0* sin 6*3 cos^ 6*3 sin6'* cos 6'* (6*0(5) + G'2('5) cos^ ^*) 
= sin6';cos6';(Go(s) + ^2(5) cos^6';) j cos 6*3 cos V^in 6^3 cos^ 6'3 

= ^sine:cose:((5o(s) +G'2(s)cos2e:). (27) 

Aj(s), i?j(s), Cj(s) and Gi are functions of Jj(s) and cos(5o(s) — ^2('5))- If 
Ai{s), Bi{s),Gi{s) and can be fitted precisely, it will not be difficult to 
determine the value of /j and 5q — 82 at definite energy. 

y4j(s), Bi{s), Gi{s) can also be written as combinations of Gkim, 
2n 

A-o = — (6G000 + 2G020 + 36*200 + ^220); (28) 

^2 = —(66*002 + 26*022 + 36*202 + 6*222), (29) 
2n 

A-i = —(66*004 + 26*024 + 36*204 + G224), (30) 
27r 

-Bo = 7^(306*000 + 106*002 + 6G004 + 156*2oo + 56*202 + 36^204)5 (31) 
15 

27r 

B2 = 77(306^020 + IOG022 + 6G024 + I5G220 + 5G222 + 36^224); (32) 
15 

4 

6*0 = 77(45Gooo + 156*002 + 9Goo4 + 156*o2o + 5Go22 + 36^024)5 (33) 
45 

4 

6*2 = 77(456^200 + 156*202 + 9G^204 + 156*220 + 56*222 + 36^224)- (34) 
45 
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The dependence of Gkim on form-factors Ii{s) and cos(5q — 5q) are very com- 
plicated. The expressions are hsted in Appendix B. 



3 Efficiency corrections 



The partial and weighted distributions in Eqs. fl24l) - (l27j) were based on a the- 
oretical integration of the partial decay rate over different angular variables. 
However, the experimental situation is slightly different from this. In general, 
the detector is not able to cover the whole solid angle and, moreover, the de- 
tection efficiency is not the same in all directions but it is a rather complicate 
function w{VL) B The partial decay rate detected in the experimental analysis 
is not that in Eq. ( |23l) but the efficiency corrected one. 



dV 



dsd cos 9* d cos 9-:! 



dr 



dct. 



dsd cos 9td cos 6'q 



(35) 



The calculation of the corrected functions corresponding to the distribu- 
tions in Eqs. flM|) - fl27|) is more tedious but it does not introduce any important 
complication. In order to ease the understanding of the procedure, we present 
a detailed calculation for d'^T /dsd cos 9* . We integrate the detected partial 
rate in Eq. fl35l) over 9^ and and we integrate separately every monomial 



cos'' 0* cos' 9^ cos™ 61* 
d^T 



dsd cos 91 



det. 



E Gkimis) jdcos9^d<Pl wiSl) cos' C cos' ^3 COS™ 9; 



k,Lm 



+ Go{s) J d COS 9 sdcf)^. w{il) cos </>* sin 6*3 cos 6'3 sin 6^* cos 6** 

+ ($2(5) J d cos 9-id(j)*^ w{Vt) COS0* sin 6^3 cos 9^ sin^* cos'^ 6^* . 

(36) 

Since the integral is on 9^ and 0* , it is possible to reexpress it in the form 
d^T_ 

det. 



E A^cos^9: + Sin 91 COS 9: 

m=0,2,4 



Ao + A2 cos^ 91 



dsd cos 9. 

where we have defined a new set of coefficients Aj, Aj given by 

= E G'Hm(s) / d COS 9-id<Pl w{p) COs' ^ COs' ^3 
fc,Z=0,2 

Aq = Gq{s) j d COS 9 ■idcj)*^ w{Vt) COS (j)*^ sin 9^ COS 9^ , 
A2 = G2{s) j d cos 9 ■id(j\ w{Vl) cos (f)*^ sin 9 ■i cos 9 ■i . 



,(37) 



(38) 

(39) 
(40) 



^Although a priori we will assume w{Q.) ~ w{d*,(f>^,93), notice that for asymmetric 
detectors the efficiency could also depend on the azimuth angle (j)3 
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It the case of perfect efficiency, w{Q) = 1, one finds Aq = A2 = and the 
different Ai become those provided in Eqs. (l28!) - (!34l) . The dependence on s 
is imphcitly assumed. Furthermore, if the efficiency depends on 6* then the 
coefficients Ai, Aj are also functions of this angle. In this case, when analyzing 
the experimental data one should compute these integrals for every 6*. The 
simplest procedure to evaluate these integrals is through the Monte Carlo 
method, where we have for instance 



rfcOsMC^(^)cOs'=0:COs'^3 - ^ E cos' C a cos' ^3, a. (41) 



In the case of a 6'*-dependent efficiency, this integral also depends on this angle 
and it must be repeated for every point in the fit analysis. 

Through a similar procedure, one also recovers the detected distributions 
corresponding to those in Eqs. fl^^ - flTTI) : 



dsd cos 03 



= + B sin ^3 cos ^3 , (42) 

det. 1=0,2 

= ^ Cfc COS V; + C COS C (43) 

det. fc=0,2 



dsd cos 

W"[^.^:]ldet. = E W^A:Cos'"^: + sin^;cos^: Wo + W^cos^Ol (44) 



m=0,2,4 



(45) 



B 



with the coefficients 

Bi = J2 GkiUs) fdcose:d<P:w{n) cos' cp^cos'^e:, 

fc=0,2 m=0,2,4 •' 

= Jdcos eid(j)l w{n) cos 0; sin e; cos ei [Go{s) + G2{s) cos^ ei] ,(46) 

= E E GkUs) fdcose:dcosesw{n)cos'e3cos'^e:, (47) 

/=0,2m=0,2,4 •' 

C = j d cos O'^d cos 6*3 (f2) sin O-^ cos O-^ sin 6** cos 6^* 

X [Go(s) + G2(s)cos2e:], (48) 
Wm = GmUs) fdcosesd<p:w{Q)cos'+'<p:cos^+'e3, (49) 

k,l=0,2 

Wo = Go{s) J d cos 93d(l)l w{n) cos^ 0; sin 63 cos^ 63 , (50) 
W2 = G2is) Jdcos 93d(l)l w{n) cos (pl sin 63003'^ 63. (51) 

As it happened before, for a general efficiency w{Q), these coefficients are not 
simply functions of the energy but they also have a residual dependence on 
the corresponding angle. 
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4 Monte Carlo Study 



Clear signals of a and k are found in BES data[T4l [151 E]- In this studies, 
measuring the tttt/ttK S-wave phase shift is tried, but because of the limited 
statistics, no meaningful results are obtained. In the J/\Ef — cuvrvr channel, 
there are resonances in the un spectrum, which will affect the S-wave phase 
in the vrvr spectrum. Its contribution to the tttt S-wave phase shift is hard to 
be estimated theoretically, which is the trouble for the measurement of the tttt 
S-wave phase shift in the J/^' — > cutttt channel. However, all these troubles do 
not exist in the — > tttt J /"^ channel, for the energy of the ttJ/'^ spectrum 
is too low and no resonances exist in the irJ/'^ spectrum. For BESII data, 
the channel \1'' rnrJ/'^, where J/^ ^ /^^/^~! is studied[6], and a global 
partial wave analysis is performed. After introducing a wide 0"^"'' background 
which strongly destructively interfere with a particle, the vrvr spectrum can be 
well fitted. The pole position measured in this channel is consistent with that 
measured in the J/\l/ — > a;vrvr channel. Though global PWA fit can obtain 
reasonable results on a particle, vrvr S-wave phase shift can not be well defined. 
The reason is that the statistics in BESII data is too low to perform a reason- 
able fit on phase shift, which is studied in this paper. In this paper, we will 
use Monte-Carlo technique to generate data with different D-wave components 
and different statistics, and then use the method proposed in this paper to fit 
the data. 

It is expected that BESIII will collect huge number of ^' events. The 
statistics of BESIII data will be about one thousand times of BESII data. For 
example, the statistics of \E'' — > vr7rJ/\E^ in a 10 MeV bin at 500 MeV in the 
vrvr spectrum in BESII data is about 1,000, if BESIII statistics is 1,000 times 
larger, we will have about 1 million statistics in a 10 MeV bin. So, in our 
Monte Carlo simulation, about one million Monte Carlo events are generated, 
and the method proposed in this paper is used to fit the data to see whether 
reasonable phase shift can be obtained or not. For the Monte Carlo data, the 
S-wave phase motion is known, so we can test the above method by comparing 
the fitted results with the input value of Monte Carlo simulation. 

In the Monte Carlo simulation, we need first to know the amount of D-wave 
component, or the percentage of D-wave component in the total Monte Carlo 
data sample. According to literature [5], the ratios of D-wave component to S- 
wave component in the 7717^,7 range from 340 MeV to 600 MeV are in the range 
from 4.7% to 31.9%, and the ratio decreases with the increase of 772,^77. In order 
to simplify the problem, we generate Monte Carlo data in the energy between 
500 MeV to 510 MeV with different D-wave component. Five independent 
Monte Carlo data samples are generated with D-component 2%, 4%, 8%, 20% 
and 45% respectively. In each data sample, the method proposed in this paper 
ia used to fit the data, then a scan on I,J=0,0 phase is done. Scan results on 
the phases of the Monte Carlo data samples are shown in figure [2l In these 
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figures, we can see that there is a minimum in the smooth scan curve, and the 
phase value at the minimum is the 1,J=0,0 phase of that case. In all these cases 
except the 45% case, the input phases are 1.17, which corresponding to 67°. In 
the 45% case the input phase is taken as 1.03 radian. The fit results are listed 
in table I. It can be seen that the fit results are quite close to Monte Carlo 
inputs. So, the phases obtained by the method of this paper are reasonable. 



Ratio of D-wave 


2% 


4% 


8% 


20% 


45% 




1.28 ±0.08 
1.52 ±0.04 


1.10 ±0.16 


1.24 ±0.08 


1.18 ±0.04 


1.04 ±0.04 



Table 1: Scan results on 6q — 6q for different D-wave components. Errors are 
only statistical. The input phase difference of Monte Carlo data are 1.03 for 
the 45% case and 1.17 for the rest. 



Statistics of data is key important in the study of the phase difference. 
In the above study, the statistics of Monte Carlo data are 600,000, 250,000, 
200,000, and 200,000 for D-wave component 2%, 4%, 8%, and 20% respec- 
tively. When the ratio of D-wave component is small, we need much higher 
statistics, otherwise likelihood function is not sensitive to the change of S- 
wave phase, or even not changed when S-wave phase is changed to any other 
value. In BESII data, there are only about 2,000 events in a 10 MeV bin when 
mTTTT is at 500 MeV. In these case, the fit is not at all sensitive to the phase 
shift difference. In the Monte Carlo study, we found similar results, that is, 
when D-wave component is 4% and statistics of the Monte Carlo data is below 
10,000, the likelihood function almost kept unchanged when we change the 
S-wave phase to other value. Therefore, the reason that we can not obtain 
a reasonable result on S-wave phase shift is that the statistics of the BESII 
data is too few. It is expected that BESIII will collect 200 times more ip' data 
in the near future. And BESIII detector has much higher selection efficiency 
than that of BESII detector. So, we will have more than 400,000 events in a 
10 MeV bin when at 500 MeV. Our Monte Carlo study shows that, if 

the ratio of D-wave component is above 2%, a reasonable S-wave phase shift 
can be obtained based on BESIII data. 



5 Conclusions 

In this paper, a method is proposed to measure tttt S-wave phase shift, and 
this method is tested by Monte Carlo data. In the Monte Carlo study, it is 
found that, even if the ratio of the D-wave component is above 2% and the 
statistics in one 10 MeV bin is above about 200,000, a reasonable results on 
TTTT S-wave phase shift can be obtained. It is expected that BESIII will collect 
enough data, so based on BESIII data, we can measure tttt S-wave phase shift 
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Figure 2: Scan results on I,J=0,0 phase, (a), (b), (c), (d) and (e) are scan 
results when D-component are 2%, 4%, 8%, 20% and 45% respectively. The 
value at minimum of the curve is the output phase. 
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in the mass region from 350 MeV to 550 MeV. 
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A Kinematics 

e' and e denote the polarizations of and J/\l/ respectively. Here we choose 



and 



M 



(IpsI, 0,0,-^3 



e± = ^(0,l,±^,0), 



ei = (0, — cos 6^3 COS0, — cos6'3 sin0, — sin 6*3), 
62 = (0, — sin 0, COS0, 0). 

P3-q = Vs{E3 + \p3\i3), 



P3 ■ ei 
p-e[ = 
q-e[ = 

P3 ■ £2 
P-e2 = 
9 ■ e'a = 
p-e+-- 
p ■ e_ : 

P ■ eo = 
q-e+-- 
g- e„ = 



= -IpsI sin6'3, 

2 I (cos 9s cos sin ^* + 7 cos 6'* sin 6'3 ) , 
^/sp'J sin 6'3, 
= 0, 

2K^|sin^; sin0, 
0, 

--V2\p:jsme:, 

2\p:^jcose:^iEs + \psr ^ 

-0, 
■0, 

1 



Ma 



q-eo = s-fdpsl + Esp)- 



I 

e'2 
e'2 



e+ = cos ^3 (cos + 2 sin 0) / -\/2, 
e_ = cos ^3 (cos — z sin 0) / v^, 
eo = --E3sin6'3/Mvi,, 
e+ = (—2 cos + sin 0) / a/2, 
e_ = (z cos + sin 0) / -\/2, 
eo = 0. 



B (^Aj/m and 

We have used the notations 
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and S = So — 62 in the definitions of Gkim and Gi [§. 

G224 = -^i2hl7P^M'' + p^i-Ie^^p,^ _^ EiV)^^ - 2hhM(P3 + EiVfiEs + P^Vh'' 



+/7^m2{m2 - (£3 + PaVf-y'^))), (52) 

G204 = — G224, (53) 
G222 = -^ip\3h'^p\P3 + EsYf-f* - Ir'^M^p'^iiM^ - SiE^ + PaVfy') 

+2l7M(3cosShM^ - ihM^ + Shp'^iPs + E^VfiE^ + PaV^^))), (54) 

G202 = — G222, (55) 

G220 = ^hpHhp' - 6cos 5/1 + 2/3), (56) 

G220 = — G200, (57) 



G024 = -j^iEa^ih + hfp^vS'^ + 2E3^{h + /6)p^v2(72(77Mp2 + v{P3(hV'^ + h+ 2/6 )7P^ + UMV)) 
-hM)-i^ + {Ir'^M'^iPs^V^'r^ - M^)p'^ + 2hMP3V'y{h{Pi^V'^'y'^ - M^)p^ + Pa'rihPs'rp'^ + IaMV))p^ 
+-P3S^(((^'3 VV - M2y2)jg2 _^ 2/6P3^VS^/5 + /6^P3S^)p* + 2lAMPaV{hV^ + h)'yp^ + M^V'^))-y'^ 
+2E3{'y{h'^P3V^'y(P3'^{V^ + 1)7^ ~ A/^)p^ + hVilrMViPa'^ (V^ + 2)7^ - Af2)p2 

+-P3S^(^6P3(3y2 + 1)7P^ + /4M(2y3 + y)))p2 ^ p^-y(ll^2yj^2 ^ 7g7^p3(2v2 + i)^M + 2/6^^3 V72)p4 

+/4My(/7M(y2 + 1) + 3/6P3V^7)p^ + /4^mV3)) - hMPaiihMV + PsihV^ + hh)p^ + hMV))^^ 

-13'^ M* + Es'^iilr^M^p-^ + 2l7MV(P3{2hV^ + h + h{V^ + 2))7p2 + /4A/y)p2 

+F2(((P32(v4 + 4y2 ^ ^-^^2 _ M^V^)l5^ + e/fif^^C^^^ + 1)7^-^5 + 6/6^P3S^)p'' + 2hMP3V{3Ie 

+l5iV^ + 2))7p2 + /42Af2y2))^4 _ 2I.j,M{ItMp^ + V(P3(/5 + h + 276)7P^ + l4MV))y^ + h'^M'^)), (58) 

G004 = ^(/6'f3VS' + h'P3''p^VS'' + 2/5/6^3 + 2/5/7 AfP3V*V^S' + 2/4/5A'/P3'p'vS' 

+2/6/7^/^3^*^75 + 2lAhMP3^p^V'y^ - I6^P3^P*"/* - h'^M'^Pa^p-^V'^'^^ + I-j'^M'^Pi^p^V'^'y^ 

+/42M2p32y2^4 _^ 2/4/7 A/2p32p2y2^4 + £;3V*V'*((7^ " l)/6^ + 2/57^/6 + /5^7^)7* " 2hh M P3^ p^^V 'y^ 
-2/5/7 A/3p3p4y^3 + 2£;33p2\/2(2/62p3y7(72 - l)p2 _ 73(75 + 7g)A/ + Is'Y^I^Mp'^ + Vil5P3iV^ + 1)7P^ 

+hMV)) + hil-rMi-r^ - l)p2 + V72(/5P3(y2 + 3)-^p2 + i^MV)))-y^ - Ir^M'^pS^ - M'^ P3^ p^V'^-y'^ 
+2E3{Ir^M^P3V{-y^ - l)p'* + hM-y(hP3^{2V^ + l){-y'^ - l)p^ + y(/5 V'(P3^(V^ + 2)7^ - A/2)p2 
+hMP3(V'^ + l)7))p2 - IsMPsiilrMV + P3{hV^ + Ieh)p^ + hMV) + P3V ■y'^ {2h'^ Ps'^ {■y'^ - l)p'' 

+/6P37(/5P3(3V'^ + 1)7P^ + 3hMV)p'^ + V(/5^V(P3^(V2 + 1)72 _ A/2)p4 + /4/5A/P3(2V2 + l)^p2 

+/4^A/V)))72 + 13^ + h'^M^p'^ + 2/3/7 A/" p2 + E3^{l3^M^ - 2/37^(/7 Mp^ + V{P3{hV'^ + h + 2h)'yp^ 
+IaMV))M + 72(/72a/2(^2 _ ^^pA ^ 2hMV^{hP3{V'^ + 2)(72 - l)p^ + -y{IsP3(2V^ + 1)7P^ + /4My))p2 
+FS^(6/6^P3^(7^ - 1)P* + /5^(P3^(V* + + 1)72 - Af2y2)p4 _^ 2/4/5 A/P3V(V2 + 2)7p2 
+6/6P37(/5P3(V2 + 1)7P^ + hMV)p^ + h^M^V^m, (59) 

G022 = ^(2/52p3V'vS« + 2/5/6P3VVS'5 + 2£3*(/5 + hf P^VS"" + 21^ P3* p^V^'^ + 2/5/6P3V'^S' 

+4/5/7 A/P3^p''V^7'5 + 2/6/7 A/P3^p''y^7'' + 4/4/5A/P3^p^y^7'' + 2IaI6M Pa^ p'^V^'r'^ + 2hhMP3^pW-i^ 
+2hIeMP3^p'^Vr' + 2Ie^P3^pS* + 2hIeP3^p^V^y^ - Sh^M^Pa'^p^V^y^ + 2/72A/2p32p4y2^4 
+2l4'^M^P3^V^y'^ + 4/4/7M2p32p2v2^4 _^ 2Iel7M P3^ p'^V^^ - Qhlj PspW'y''^ + 2/4/6 A/P3^p2y7^ 
+2E3^{h + /6)pV2(/7Af((V2 + 1)^2 ^ ^^p2 _ ^j^j^ ^ ^^^^^^^^^^2 ^ 3^^2 ^ -^^^2 ^ 27(/5P3(F^ + 1)7P^ 
+/4A/y)))73 - 3/72A/4pS^ + 2£3(7{(/72p3\/((V2 + 1)72 + l)A//2 + /77(/6{72y4 + (5^2 _^ ^)y2 ^ 2)P32 
+/5V2(P32(3(\/2 + 1)^2 ^ _ 3A72));y^ _^ pgy^2(^g2((3^2 _^ ^)^2 _^ 3)p^2 ^ J^J^(^2'y^V'^ + (67^ + 2)^2 + ^2 

+1)P3^ + /5^v2(2P32(y2 + 1)^2 _ 3M2)))p4 + UM P3V {h M {{'iV'^ + 1)7^ + 1) + 2P3V7(/5(2V2 + 1)72 
+/6((y2 + 2)72 + l)))p2 + 2l4'^M'^P3V^y^) - l3MP3{hP3{iV^ + 1)7^ + 1)P^ + 2Vy{{l7M + hP3V'y)p^ 
+/4A/)))7 - 2/3^Af* + 2/3/7 A/'' p2 + 6cos5A/(-/2P3^V((/7M\/ + P3(/5V^^ + /6)7)p^ + /4My)7* 
+-B3^(/5 + /6)pV2(/i - /2yS^)7^ + E3P3ihi{l7MV + PsihV^ + h)'y)p^ + I4MV) - hViy'^ihMiV'^ 

+ I)p2 + y(P3(2/5y2 + 75 + 3/6)7P^ + 2I4MV)) - l3M))-f^ + hhM^ + £;32(/i(72(/7A/p2 

+V(P3(/5V^^ + /5 + 2/6)7P^ + hMV)) - I3M) - /2y27^(7^(/7A/p2 + y(P3(3/6 + h(V^ + 2))7p2 + UMV)) 

-I3M))) + E3^(2l3^M^ - 2/3(/7A/((y2 + 1)^2 ^ -^)^2 ^ y^^j^p^^i^y2 ^ 3)^2 _^ ^)^2 ^ 2-f{hP3{V^ + l)7p' 

+/4A/V)))A// + 72(2/72A/2(y2^2 _^ ^)^4 _^ 2/7 A/V^ (P37(/6 (5^^^ + 7^ + 3) + /5(7^V''' + 47V2 + y2 

■^The Mathematica notebook and fortran programs can be obtained from xi- 
aoly@pku.edu.cn and cillero@ifae.es. 
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+7^ + + hMViiV'^ + 1)7^ + l))p2 + y27(7((2P32(y4 + 4^2 ^ -^^^2 _ 3j^2y2)j^2 _^ 2/6^3^(7^^* 

+ (77^ + ^ 4^2 _^ 2)75 + Gh^Ps'aV^ + 1)7^ + 1))P'' + 2lAMP-iV{2h{V'^ + 2)7^ 

+-r6((V^ + 5)7^ + l))p2 + 2/4^mVS)))), (60) 

G002 = -^{2I^P-i^p^VS'' + 2/5/6^3^*^^6 + 2/62P3VV7' + 2/5/6P3V*V276 +4/5/7MP33p4y375 

+2/6/7MP3^p''V^7'5 + 4/4/5MP3^P^V^7^ + 2/4/6MP33p2V^^7^ + 2/6/7MP3^p''y7^ + 2/4/6 MP3^p2y7^ 
-h^P-i'^pS'- + 2/5/6P3V*V274 - 3/52M2p32p4v2^4 + 2/7^ A/2p32p* V^* + 2/42M2P3V274 
+4/4/7M^P3VV27'' + -B3*P''V*((272 - 3)/62 + 4/57^/6 + 2/5^^)7'' - 4/6/7 A/P32p'*y^3 _ 6/5/7M2P3p*y73 
+2hhMPi^p'^V-i^ + 2E-/p^V'^{2h'^P3p^V{V'^ + 1)7=* - 2/3(/5 + h)M + h{l7M({V^ + 1)7^ - 2)p2 
+Vi{hPi({V'^ + 3)7^ - 5)p2 + 2hMV^)) + h{l7M{{V^ + l)-f^ + l)p^ + V-yihPaii^V^ + 5)7^ + l)p^ 
+2/4MV7)))73 - 3/7^M4p''72 - 3/7^M2p32p4v2^2 ^ 2P3(7({/72p3y((F^ + 1)7^ - 2)A/2 + /77(^6(7V^ 

+5(7^ - _ -^)pg2 _^ j^^2(pg2(3(^2 _^ ^)^2 _^ ^) _ 3A^2))j^^ ^ pgy^2(^g2((3y2 ^ -^)^2 _ 3)pg2 

+/5/6(27^V'' + (57^ + 2)y2 + -y2 + i)P32 + jg2y2(2P32(y2 _^ ^^^2 _ 3M'^)))p'^ + UMP^VihMiiW'^ + 1)72 

+ 1) + 2P3y7(/5(2y2 + 1)^2 ^ j^^(y2 ^ 2)-y2 _^ ^)))^2 ^ 2/42M2p3\/S^) " /s A/P3 (/6P3 ((F^ + 1)7^ + l)p^ 

+2V'y{{hM + hP3V-i)p'^ + l4M))h + 2h'^M^ + ilr^M^p'^ + HhhM^p^ 

-GcosSMihPa^ViihMV + PsihV^ + hh)p'^ + hMV)'y'^ + Ea^ih + h)p^V^il2V^'y^ - h)y^ 

+E3P3(l2V{'y^(l7M(V^ + l)p2 + V{P3(2hV'^ + h + Shhp^ + 2hMV)) - I3M) 

-hiihMV + PsihV^ + hh)p^ + l4MV))y^ + /iA/='{/7p2 + /g) + Es^ {hV^ {■^'^ {h M p^ + y(P3{3/6 

+h{V'^ + 2))7p2 + /4MV)) - /3A/)72 + h(hM - y\l7Mp^ + ViPaihV^ + h + 2Iehp^ + hMV))))) 

+E3^(2h'^M'^ - 2/3(/7M((y2 + 1)^2 ^ ^^p2 _^ y^(/gP3((v2 + 3)72 + l)p2 

+27{/5P3{v2 + l)7p2 + UMV)))M + 72{/72a/2(2V^72 - l)p4 

+2/7My(P37{/5{7^V* + {472 + l)y2 + ^2 _^ ^) ^ /6{(572 - 3)y2 + ^2 _ 3))^2 

+IaMV{{V'^ + 1)7^ + l))p2 + yS{7({2P3^{V^ + 4^2 + 1)^2 _ 3j^,^2^2)^^2 

+2/6P32(7^V'' + (77^ + 1)^2+472/5 + 6/62p32{(y2 + 1)^2 _ 2))p4 

+2/4MP3y(2/5(y2 + 2)72 + /6((y2 + 5)72 + l))p2 + 2/4^4/2^2^)))))^ (gl) 
G02O = -^(^52P3V*VS' + /62P3V*FS' + 2/5/6P3V*VS' + i?3*(/5 + h? P^VS^ + 2I5I7 M Pa^ p'^V^^'' 

+2hl7MP3^p'^V^-^^ - Qcos ShUMP-i^p^V^"/^ + 2/4/5A/P3^p2v^3^5 _ QcosShhMPs^ p^V^^^ 

+2l4hMPz^p^V^'y^ + 2/62p3''p*yS* + 2/5/6P3''p''v274 + /72M2p32p4 1/2^4 _^ g^^2^2pg2y2^4 
+/42M2p32y2^4 _ Qcos ShhM'^Pi^V'^-l'' - 6cos5/2/7M2p32p2y2^4 _,_ 2/4/7M2p32p2l/2^4 
+2/6/7 A/P3V*y7^ - 6cos5/2/6MP33p^V7^ + 2/4/6MP3 ^^y^S _^ 2£:3^{/5 + /6)p2 V2{/5P3p2v3^3 

+/6P3p2 + /5P3p2 V73 + /6P3p2y73 + /7MPV272 + /4MyS^ + hPsp'^V-y + /7A/p2 

-/3M + 3cos5M(/i - /2V272))^3 _^ ^g2pg4^4^2 ^ 2E3P3{h'^ Pa^ p'^V^^ + h'^ Ps'^ P^V^^ 
+2hIeP3^p'^V^-y'^ + /52P32pVS^ + /6^P3^P''V^S^ + 2/5/6P3 + 2/5/7MP3p4yS* 

+2/6/7MP3p''y*7'' + 2/4/5MP3p2y474 + 2/4/6MP3p2y47'' + /5/7MP3p''yS'' + /6/7MP3p*yS* 
+/4/5MP3PV27'' + /4/6MP3p2y2^4 ^ 772^2^4 y3^3 ^ 2/62p32p4y3^3 ^ 2/5/6P32p*yS^ 
+9/2^71/2^3^3 _,_ /42j^^2^3^3 _,_ 2UI7 p'^V^-^'^ + /62p32p''y73 + /5/6P32p''y7^ + hh M P^pW^ l'^ 
+3Iel7MP3p*V^y^ + 2/4/6MP3p2y2^2 ^ 7^2j^^2^4y^ ^ /62p32p*y7 - dhhM^V^ 
+/4/7M2p2y-y + If^I^MPip'^ - /3M(/6P3(y272 + 1)^2 _^ + hP3V'~l)p^ + /4M)) 

+3cos<5A/(/i(/6P3(y272 + i)p2 ^ V^dhM + /5P3V7)p2 + /4M)) - /2y7(^7A/{2y2^2 _^ ^)^2 

-/3A/ + y7(/6P3{(2y2 _^ ^)^2 _^ 2)p2 + 7(/5P3(2y2 _^ ^)^^2 _^ 2/4My)))))7 - g/i^M* 

_732m4 - /72A/V + 2/3/7MV + 6cos5/iM4(/3 - /7p2) + £;32(/52p32p4y67^ + Ps^ p'^V'-y'' 

or T d2 4t/6 6 , «7-2o2 4t/4 6 , «7-2o2 4,^4 6 , ot t d2 4t/4 6 , t2o2 4t/2 6 
+2/5/6P3 P V 7 + 4/5 P3 p y 7 + 4/6 P3 p y 7 +8/5/6P3 PV'7 +/5P3PV^7 

+/62P32p''y2^6 _^ 2/5/6P32p''y2^6 + 2/5/7 A/P3P''VS'^ + 2hl7 M Pip'^V'^'^ 
-6cos .5/2/5 A/P3p2y575 + 2/4/5MP3p2y5^5 _ 6cos 5/2/6MP3p2y575 + 2/4/6MP3p2y5^5 
+4/5/7 A/P3p''y3^5 _,_4j-g^^jy^pgp4y3^5 - 12cos5/2/5A'/P3p2y^7^ +4/4/5MP3p2y3^5 
-12cos5/2/6A/P3p2y3^5 ^ 4/47gMP3p2y3^5 ^ /72^^2^4y4^4 ^ 2/62p32p4y474 
+2/5/6P32p''y''7'* + 9/22M2y''7* + /42M2y*7-* - 60085/2/4^/2^47-* - 6cos <5/2/7M2p2y4^4 
+2/4/7 A/2p2y474 + 4/62p32p4y2^4 _^ 4^^ jgPg2^4y2^4 _^ 2/5/7MP3p4y3^3 _^ 4/g/7MP3p4y373 
-6cos5/2/6MP3p2y^7^ + 2/4/6MP3p2y3-y3 + 2I5I7M P3p'^V'r^ + 2IqI7 M Pip-^V^i^ 
+2/72M2p4y2^2 _|_ ^g2pg2^4^2^2 _ g^^^ gj^j^j^2^2y2^2 _^ 2l4l7M^ p'^V^y^ + 2/6/7MP3p'*y7 
+/72M2p4 + 9/i2a/2 + /32m2 - 2/3 A/(/7A/{y272 + l)p2 + y7(/6P3((y2 + 1)72 + l)p2 
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+7(^5f'3(^ + 1)7P - 3cos<572MF + I4MV))) + 6/1^(008 S{l7M(V^Y + l)p - hM 

+V'rihP3iiv'' + 1)7^ + + -rihPsiv^ + ihp^ + UMV))) - shMV-/'^))), (62) 

-6cos ShhMPz'^p^sV'-^-i^ + 2/4/5 MP3^p2gy3^5 - 6cos ShhMP^^ sV^-^^ + 2/4/6A/P3^p2^y37^ 

+2/6^P3''p*sy^7'' + 2/5/6P3''p''sy^7* + 9/2^M2p32sy2^4 _^ j-^2^2p^2^y2^4 

-6cos5/2/4A/2p3^sF27'' - 6cos .5/2/6 A/Ps^p^^y^a + 2/4/5 A/P3^p^sV^7^ 

+2E3^{h + h)p^V^i-hMs + 3cos5M(/i - /ayS^)^ + V^7(P3(/5s(y2 + 1)7^ + /6s((V^ + 1)7^ + l))p^ 

+/4A/sy7))7^ + h^Ps'^p'^s-y^ + 2£;3P3s(P3^V7(^5^V2(V2 + 1)7'' + /5/6(27^ + 2(7^ + + 1)^2 

+/6^((V* + y2)^4 ^ (2^2 + 1)^2 _^ ^))^4 _^ MP3(/4V'2(/5(2y2 + 1)7^ + /6(2yS^ + 7^ + 2))7^ 

-hihV^l' + /eVS^ + -fe) + 3cos5(/i(/5yS^ + ^e^S^ + h) - l2V'^-i'^{h{2V^ + 1)7^ 

+h{2V'^-1^ + 7^ + 2))))p2 + M2y7(9/2 + li^V^'l'^ - hh + /i(3cos SI4 - 9/2) 

+3cos<5/2(/3 - 2l4,V'^'y^)))'y + dh'^ M'^s + M'^s - 6cos 5/i/3M''s + /72M2p-l5(4M2 

+ (S3V^2^2 _^ pgy^2 _^ ^^^2) ^ 2hMp^{-l3M{Ezs{E-iV'^'y'^ + PaV-y'^ + E3) - 2M'^s) 

+ (P3 + EsVhiEsV'^j^ + PsV-y'^ + B3)(-f6s(P3 + ^3^7^ + P3)p2 + sVy{h{E3 + P^Vhp^ + I4M)) 

-3cos5M(hsViP3 + EaVh'^lEaV'^j'^ + P^V-i^ + E3) - I^_{E3s{EzY'^ y^ + P-iVy^ + B3) - 2M2s))) 

+E3^(9/i^sM2 + /32sm2 - 2/3V7(P3(/5^s(v2 + 1)7^ + /6S({V2 + 1)7^ + l))p^ 

+ (/4 - 3cos<5/2)MsV7)A/ - Ux{3l2M sV'^ + cos5s(/3M - V7{/6P3((V^^ + 1)7^ + l)p^ 

+7(^5^3(^2 + i)7P^ + UMv))y)M + yS^(P3^s(/52(y* + + 1)^4 

+2/5/6(^2 + (y4 + + 1)^2 _^ 2)72 + /g2 ((y4 _^ ^ -^)^4 ^ 2(y2 + 2)72 + i))p-i 

-2(3cos 5/2 - /4)MP3y7{/5s(y2 + 2)72 + /6s((y2 + 2)72 + l))p2 

+ (9/22 - 6cos 5/4/2 + /4^)M2sy2^2)-))^ (g3) 
Go = -^(2p27{-2/72p2M'* - /3/7A/4 _ 2/5/7 P3p2y^M3 - hhPiVyM^ + /4/7P3V272M2 

+/72P32p2y2^2^2 _^ I^I^p^3 p'^V^-f^ M + IqI-jPs^ p^V^-y^ M + hhPs^Vy^M + hh P^^ p'^Vy'^ M 

+hl7P3^p'^VyM + ScosSihV^-yih - l2V^-f^)E3^ + ihilrM + 2IeP3V'y) - hV^y^i^M + 3l6P3Vy))E3^ 
+ {h{hyP3^ + I7MVP3 + hM'^V^y) - /2P3V72(/7A/(y2 + i) + 3IeP3Vy))E3 
-hPs^V-r^ihMV + /6P37) + /iA/2(/7A/ + /5P3V7))A/ + E3*h{h + h)p^V^y* + le^ P3'^ P^V^y* 

+hhP3'^P^V^1* + /6^P3*pS^ + £;3='v27{/5/7A/72p2 + h'^P3Vy{{V^ + 3)72 + l)p2 - /g/gM 
+/6(/7Af{(y2 + 1)^,2 + l)p2 _^ y^2(^^p^(y2 ^ ^-^^^2 ^ l^MV))) + E3^ {I-j'^ [V'^ y'^ + l)p2 

-hMihM + 2hP3Vy) + 3hP3V'^y'^{hP3{{V^ + 1)7^ + l)p^ + 7(^5P3(V2 + l)7p2 + I4MV)) 

+/7A/y7(/6P3((5V2 + 1)72 + 3)p2 + 7(/5P3(2V2 + l)7p2 + I4MV))) + £;3(/7^A/2p3y((y2 ^ ^^^2 ^ -^^^2 

-hMile-yPa^ + /7A/VP3 + /5A/2y2^) + /6P3 V72(/6P3((3y2 + 1)^2 ^ 3)^2 
+7(4^3(3^2 _^ ^)^^2 _^ 37^A,/y)) + IjMy{hP3^{y'^V^ + (572 + 1)^2 _^ 2)p2 

+ y(/5y(P3^(V^^ + 2)72 - 2A/2)p2 + /4MP3(y2 + 1)7))))), (64) 
G2 = -jgj(2p27(-/72p2M4 - /5/7P3PV7A/3 _ /3(/7M + /5P3y7)A^^ + /4/7P3^ ^2727^2 

+/72p32p2v2T,2M2 + hl^P-/ p^V^y^ M + hhP3^Vy^M + 2hI-jP3'^p^Vy^M + E3'^h{h + h)p^V*l'^ 
+h^P3^ph* + hhPs^p'^V^'y* + E3\-y\lr^M'^p' + ShPsV^jiPaihV^ +h + 2h)yp^ + hMV) 
+/7A/y(P3(2/5V2 + /5 + 2/6(y2 + 2))7p2 + UMV)) - hMiljM + 2/6P3V7)) 
+B3V27(72(/5(/7M + /6P3l^(l^^ + 3)7)p2 + h{2lT M p^ + V{ihP3yp^ + I4MV))) - hhM) 
+E3{-l{hV(hP3^{3V'^ + 1)7^ + I7MV{P3^{V^ + 2)72 - A/2))p2 + P3^(2(/72l/A/2 +/677P3(2V2 + l)yM 
+2/62p32y72)p2 + UMVihMiV'^ + 1) + 3/6P3 V7))) - l3M{hyP3^ + /7A/VP3 + /5A/2y2^)))). (65) 
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